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Abstract 

The renormaUzation group method is apphed to the study of homogeneous and flat Friedmann- 
Robertson- Walker type Universes, filled with a causal bulk viscous cosmological fluid. The starting 
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^J^' point of the study is the consideration of the scaling properties of the gravitational field equations, 

of the causal evolution equation of the bulk viscous pressure and of the equations of state. The 
^ ' requirement of scale invariance imposes strong constraints on the temporal evolution of the bulk 

viscosity coefficient, temperature and relaxation time, thus leading to the possibility of obtain- 
ing the bulk viscosity coefficient-energy density dependence. For a cosmological model with bulk 
viscosity coefficient proportional to the Hubble parameter, we perform the analysis of the renor- 
malization group flow around the scale invariant fixed point, therefore obtaining the long time 
behavior of the scale factor. 
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I. INTRODUCTION 



The renormalization group (RG) has been proved to be a powerful tool in statistical 
physics, quantum field theory, critical phenomena and chaotic dynamic ^. The essence 
of the renormalization group method is to extract structurally stable features of a system, 
which are insensitive to details. Many physical systems in field theories, polymers etc. 
exhibit universal scaling functions and critical exponents in the limit A/C,, where A is some 
ultraviolet cut-off and ^ is a temperature-dependent correlation length. The RG group 
method is the main tool with which to elucidate this universal behavior. On the other hand, 
from a purely mathematical point of view, the RG is regarded as a means of asymptotic 
analysis. Long-time asymptotics of non-linear partial differential equations using RG ideas 
have been studied recently by Bricmont, Kupiainen and Lin |Q and by Chen, Goldenfeld and 
Oono P|. In this formulation, the RG transformation is the integration of the equation up to 
a finite time, followed by a rescaling of the dependent and independent variables. Therefore 
the problem at infinite time is reduced to the problem at finite time. A fixed point of the 
RG transformation corresponds to a scale-invariant solution of the differential equation, and 
the long-time behavior of the solution can be obtained by studying the flows around flxed 
points. The RG method for reduction of evolution equations has been formulated in terms 
of the notion of invariant manifolds in [^], and it has been shown that the perturbative 
RG method constructs invariant manifolds successively, as the initial value of the evolution 
equation. 

The RG method has been also intensively applied in the study of cosmological models 
and of the critical behavior of the gravitational collapse of a radiation fluid. Koike, Hara 
and Adachi presented a scenario based on the renormalization group that can explain the 
universality and scaling observed in a numerical study of gravitational collapse of radiation 
fluid. This approach for understanding and analyzing critical behavior in the gravitational 
collapse has been generalized, to a general framework, for perfect fluids with pressure pro- 
portional to density and with arbitrary 7 ||^. Critical phenomena occur for 1 < 7 < 1.88 
and the uniqueness of the relevant mode around a flxed point is estabhshed by Liapunov 
analysis. 

The averaging hypothesis tacitly assumed in standard cosmology was discussed by using 
RG methods by Carfora and Piotrkowska 0. A critical behavior, which can be related to 
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the formation of sheet-like structures in the Universe, is obtained, and exphcit expressions 
for the normahzed Hubble constant and for the scale dependence of matter distribution 
are also provided. The observable scaling properties in Newtonian cosmology have been 
considered in [^. The quantitative scaling properties have been analyzed by calculating the 
characteristic exponents around each fixed point, thus leading to the possibility of a fractal 
structure of the Universe beyond a horizon scale. The RG method has been applied to the 
Einstein equations in cosmology by Iguchi, Hosoya and Koike who carried out a detailed 
analysis of renormalization group flows in the vicinity of the scale invariant fixed point in the 
spherically symmetric and inhomogeneous dust filled Universe. The first order solution of 
long- wavelength expansion of the Einstein equations has been improved in [|l^] . By assuming 
that the RG transformations have the property of Lie group, the secular divergence caused by 
the spatial divergence terms can be regularized and absorbed to the background seed metric. 
The renormalized metric describes qualitatively quite well the behavior of the gravitational 
collapse in the expanding Universe. The back reaction effect due to the inhomogeneities 
has been incorporated into the framework of the cosmological perturbations, using the RG 
method, by Nambu [0. The RG equations, which determines the dynamics of the space- 
time, have been derived in a gauge invariant manner. 

The basic techniques and ideas for the construction of RG theories for the study of 
non-equilibrium processes in inflationary and semiclassical cosmology and stochastic gravity 



have been recently reviewed in [1^]. The complete study of the flow of the renormalization 
group in the Einstein-Hilbert truncation of quantum gravity in four dimensions and the 
classification of the RG trajectories have been done in [|l^]. Cosmologies with a time- 
dependent Newton constant and a cosmological constant have been investigated, from the 
RG point of view, by Bonnano and Renter |]16|. The scale dependence of G and A is governed 
by a set of RG equations, which are coupled to the Einstein equations in a consistent 
way. This model could explain the data from recent observations of red-shift of Type la 
Supernovae without introducing a quintessence field. The cosmological implications of the 
running of the cosmological constant and whether the scaling dependences of G and A spoil 
primordial nucleosynthesis have been discussed, using RG techniques, in |T7[] . 

It is the purpose of the present paper to apply the RG group techniques to the analysis 
of the bulk viscous cosmological models. 

Dissipative processes are supposed to play a fundamental role in the evolution and dy- 
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namics of the early Universe. Over thirty years ago Misner |T8[ suggested that the observed 
large scale isotropy of the Universe is due to the action of the neutrino viscosity, which was 
effective when the Universe was about one second old. Bulk viscosity may arise in different 
contexts during the evolution of the early Universe. Some physical processes involving vis- 
cous effects are the evolution of cosmic strings, the classical description of the (quantum) 
particle production, interaction between matter and radiation, quark and gluon plasma, 
interaction between different components of dark matter etc. ||19| . 

The first attempts at creating a theory of relativistic dissipative fluids were those of 
Eckart ||2^ and Landau and Lifshitz |2ll]. These theories are now known to be pathological 
in several respects. Regardless of the choice of equation of state, all equilibrium states in 
these theories are unstable and in addition signals may be propagated through the fluid at 



velocities exceeding the speed of light |22]. The problems arise due to the first order nature 
of the theory, i.e., it considers only first-order deviations from the equilibrium leading to 
parabolic differential equations, and hence to infinite speeds of propagation for heat flow and 
viscosity, in contradiction with the principle of causality. Infinite propagation speeds already 
constitute a difficulty at the classical level, since one does not expect thermal disturbances to 
be carried faster than some (suitably defined) mean molecular speed. Conventional theory 
is thus applicable only to phenomena which are quasi-stationary, i.e., slowly varying on 
space and time scales characterized by mean free path and mean collision time This 
is inadequate for many situations in high-energy astrophysics and relativistic cosmology 
involving steep gradients or rapid variations. The deficiencies can be traced to the fact that 
the conventional theories (both classical and relativistic) make overly restrictive hypothesis 
concerning the relation between the fluxes and densities of entropy, energy and particle 
number. 



A relativistic second-order theory was found by Israel [22| and developed by Israel and 



Stewart |23[, into what is called transient or extended irreversible thermodynamics. In 
the second order theory the deviations from equilibrium (bulk stress, heat flow and shear 
stress) are treated as independent dynamical variables leading to a total of 14 dynamical 



fluid variables to be determined. However, Hiscock and Lindblom and Hiscock and 
Salmonson |]25| have shown that most versions of the causal second-order theories (the so- 
called truncated models) omit certain divergence terms. When these terms are added to the 
evolution equation of the bulk viscous pressure, one obtain the full causal theory, with well 
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behaved solutions for all times. Causal bulk viscous thermodynamics has been extensively 
used for describing the dynamics and evolution of the early Universe and some astrophysical 
process ^|. 

The fact that the Universe has some hierarchical structure, and that the scale invariant 
solutions of the gravitational field equations plays an important role is well known |]TT|. 
The homogeneous Einstein equations are scale invariant. Using this property, we apply the 
RG method to study the long-time asymptotics of the general relativistic field equations, 
describing the evolution of a homogeneous causal bulk viscous fluid filled Universe. For the 
RG method we closely follow the approach presented in fill . 

The renormalization group approach shows the particular role played by the exponent 
s of the density in the equation of state of the bulk viscous coefficient. The gravitational 
field equations are invariant under scale transformations only for the value s = 1/2 of the 
exponent. Hence this method provides an effective and efficient tool for investigating general 
properties of cosmological models. 

The present paper is organized as follows. In Section II we give a brief and formal 
presentation of the renormalization group method and its applications to the study of the 
asymptotics of differential equations. The scaling properties of the basic equations of the 
causal thermodynamics in a homogeneous space-times are investigated in Section III, and 
their implications on the equation of state of the cosmological causal bulk viscous fluid 
are considered. The flow in the RG around the flxed points of the renormalization group 
equations is analyzed in Section IV. Finally, in Section V we discuss and conclude our results. 



II. RENORMALIZATION GROUP AND DIFFERENTIAL EQUATIONS 

In this Section we briefly review, from a formal and abstract mathematical point of view, 
the formalism of renormalization group, as applied to the study of systems of ordinary or 
partial differential equations. The presentation closely follows 0| and |p, and for further 
details and speciflc applications the reader is referred to these papers. 

The basic idea in the RG formalism is to consider the time evolution of initial data as 
a flow of a renormalization group transformation. Suppose we are interseted in the time 
evolution of n unknown functions y = {yi,y2, ■■■,yn), which are real- valued functions of the 
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time t and of a spatial coordinate x, and satisfy a system of partial differential equations 



dt dx' 

The time evolution of the system (|l|) is determined once one specifies the values y {t,-) 
at the initial time t = 1. The space of functions y {t, ■) at a fixed time is called the phase 
space r of the system, and assume it is a real vector space. As a first step in introducing 
the renormalization group transformation, we define the scaling transformation S {L,d,f3), 
which depends on the real parameters L E R,d = (ai, 02, a„) E R, (3 E R 

S {s, a, (3) : (t, x) ^ # (t, x) ^ e^'^y, (e-^t, e"^^x) . (2) 

If ?/ is a solution of the Eq. (|I]), and if with a suitable choice of the constants a, (3 and for 
an arbitrary scale parameter L (close to zero), the new function x) is also a solution 

of Eq. (|I|), then the system (|I]) is said to be invariant under the scaling transformation (^. 
In the following we denote, for simplicity, y^^'°''^'^ by ^y], and assume that generally L > 0. 

The renormalization group transformation 3?^ is defined as a transformation acting on 
the phase space V of the system (|lD [§]: 

: (0 = (1, (0 = W (1, = e-'^y. (e-^ 6"^^^) • (3) 

More exactly, Y is given by developing the initial data y = Y (x) from t = 1 to 
t = by Eq. (||), and then rescaling the spatial coordinate x and the unknown functions 
yi. Because of the invariance under scaling, forms a semi-group with the property 
^Li+L2 = 3^L2 ° ^Li- From this property one can find the infinitesimal generator 3? of 3?l, 
given by the derivative of with respect to L at L = 0, = lim^^^o ^"^f ^ , leading to 

= exp (^m^ . 

A fixed point of the renormalization group I L G R} is a point in F satisfying the 
equation 3?^ (Y*) = Y* for all L > 0. Equivalently, it can also be defined as 3? (Y*) = 0. A 
function y is called self-similar with parameters (a, /5) if it satisfies y(t,x) = y^^'°''^\t,x). 
Each self-similar solution ygg of Eq. ([^) with parameters (a, (3) is related to a fixed point Y* 
of 3fti by yss {t,x) = [x (-t)"^] §. 

As a next step we introduce the tangent map of 3?^ at Y. It is defined by T^^ = 
lim^^o . For a fixed L and Y, T^y is a linear operator, related to the renor- 

malization group by the following formal relation: 3^^ {Y + K) = ^l+Tl,y {K)+0 {KY + .... 



At a fixed point Y* of 3?^, the tangent map is Ti,y* = Tl (K) = hm,^o ^^(^*+^-^)-^^(^*) . 
The tangent map at the fixed point forms a semi-group. In terms of the infinitesimal gen- 
erator T = lim^^^o the tangent map can be represented as = exp (^LT^ 

The flow of the renormalization group transformation near the fixed point Y* is deter- 
mined by the eigenmodes of T, with eigenvalues /e G M, C, defined as Y*K = kK. A mode 
with Re A; > is called a relevant mode, and it is tangent to a flow diverging from Y*. A 
mode with Re A; < is tangent to a flow converging to the fixed point. Modes with Kek = 
are called marginal. If there is a finite number of relevant modes, or no relevant mode, then 
the long-time behavior of nonlinear partial equations can be predicted from the study of the 
perturbations around the fixed point. 

To make the abstract approach presented above more clearer from the point of view of 
the physical applications, we discuss the RG method for an equation of the type |Q 

y = y" + F {y, y', y") , (4) 

where a dot and a prime denote the derivatives with respect to the variable t and x, re- 
spectively. We are interested in the asymptotics of the solution of the form y {t, x) ~ 
as t — > oo. The RG methods transforms the problem of large time limit 
into an iteration of a fixed time problem, followed by a scaling transformation. The scale- 
invariant solution emerges then as a fixed point of a map in the space of initial data, the 
RG map, and the stability analysis becomes the analysis of the stability of the fixed point 
under the RG 1. 

Let us fix for the Eq. some (Banach) space of initial data S. Next, we chose a number 
L > 1 and set V (t, x) = L'^y {L'^t, Lx), where y solves Eq. (|) with the initial data / G E. 
The RG map 3? : S — > S is in this case (3?/) (x) = V (1; V {t, x) satisfies the equation 

d^y\t,x) _ d'^y\t,x) , W/^W,, , df{t,x) d'f {t,x) \ 
~^r~ - dx^ + ^ I ^ dx ' dx^ I' 

where F {a,b,c) = L'^^^'F {L '^a,L ^ '^b,L ^ "c). The renormalization group transforma- 
tion depends on L and the functional form of F. The semi-group property of the trans- 
formation implies ^l",f = 3^l,F"-i ° ••• ° ^l,Fl ° ^l,f- Letting t = L^", the long time 
asymptotics of the solution of Eq. (^) can be obtained by iterating the fixed time problem: 
y {t,x) = {^Lu^pf) (t-V2^), for t oo. 



If there is an a so that we can find two functions F* and /* so that F^n F* and 
^L",Ff — ^ /*, so that F* satisfies the scale invariant equation y = y" + F* and /* is the 
fixed point corresponding to this equation, defined as ^L,F*f* = /*, then the asymptotics 
of the initial problem are given by y (t,t^^^^x^ ~ t~°'/'^ f* [x). Therefore the long time 
behavior of the differential equation is determined by the fixed point in the RG group, and 
the stability analysis can be also performed by studying the perturbations around this point. 



III. RENORMALIZATION GROUP TRANSFORMATIONS AND THE EQUA- 
TIONS OF STATE 

The energy momentum tensor of a bulk viscous cosmological fiuid filling a fiat FRW 
Universe, with line element: 



ds'^ = dt^ - a^{t) [dx^ + dy"^ + dz^) 



(6) 



is given by [^2|, 



T^ = {p + P + Ti) n,M'-(p + n)5i 



k 
i 1 



(7) 



where p is the energy density, p the thermodynamic pressure, 11 the bulk viscous pressure 
and Ui the four velocity satisfying the condition WjU* = 1. We use units so that SttG = c = 1. 
The gravitational field equations, together with the continuity equation Tj^f, = 0, imply 



3H^= p, 
2H + = - n, 
p + 3{p + p)H= -3HU, 

where H = d/a is the Hubble parameter. 

The causal evolution equation for the bulk viscous pressure is given by [^, ^ 



(8) 
(9) 
(10) 



rtl + U = -3^H-KuhH+^-^-^ 



(11) 



where T is the temperature, ^ the bulk viscosity coefficient and r the relaxation time. Eq. 
(|lT|) arises as the simplest way (linear in 11) to satisfy the H theorem ( i.e., for the entropy 
production to be non-negative, 5"*^ = H^/^T > 0, where 5* = eN^ — -^m* is the entropy 
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flow vector, e is the entropy per particle and A^* = nw* is the particle flow vector) p3| , p4 . 



When the the condition 



dt I 



<< 1 holds, the additional terms to the evolution 
equation are negligible in comparison with 3C,H, and the full equation leads to a truncated 
equation with reduced relaxation time and reduced bulk viscosity p3| . 

We also suppose that the bulk viscous cosmological fluid obeys the barotropic equation 
of state 

p = (7 - l)p, 1 < 7 < 2. (12) 

With this choice the conservation equation (0) can be formally integrated to give 

p{t) = ^{po + J a'^^'aUd?! , (13) 

where po ^ is a constant of integration. 

We consider the following scale transformations: 

t -> Lt, (14) 

a{t) ^ a\t) = L"a(Lt), (15) 

where L > 1 is the parameter of the scale transformation and n is a constant. With the use 
of Eqs. (p!^-(p^), from the Hubble parameter definition and from the field equation (^, it 
follows that the Hubble parameter and the energy density are scaled in the following way: 

H{t) H{t) = LH{Lt), (16) 
pit)^^p\t)=L'piLt). (17) 

From equations (||) and (|10]) we can now easily find the scaling properties of the bulk 
viscous pressure and of the scale factor respectively: 

H(t) ^ H (t) = L^H (Lt) , (18) 

(L) 2 

a{t) a {t) = L~^a{Lt). (19) 

Eq.(|TT]) is scale invariant with respect to transformations (P^)-(P^), if the relaxation time 
r, the bulk viscosity coefficient ^ and the temperature T obey the following scaling laws: 

r{t) ^^r\t) =L-\{Lt), (20) 

at)^ ut) = LaLt), (21) 

T{t) ^ T{t)= L^T {Lt) . (22) 
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where m is a constant. The causal evolution equation does not impose a specific scaling law 
for the temperature. 

In order to close the system of field equations (||)-(|11|) we have not only to give the 
equation of state for p, but also specify T, r and C,- There are many proposals in the 
physical literature for equations of state for the relaxation time, temperature and bulk 
viscosity coefficient 1 , [0, ^ . 



Murphy |^6|], Belinskii and Khalatnikov and Belinskii, Nikomarov and Khalatnikov 



^ have analyzed the relativistic kinetic equation for some simple cases. Their results show 
that in the asymptotic regions of small and large values of density the viscosity coefficients 
can be approximated by power functions of the energy density, with definite requirements 
on the exponents of these functions. For small values of the energy density it is reasonable 
to consider large exponents, equal, in the extreme case, to one. For large p, the power of 
the bulk viscosity coefficient should be considered smaller (or equal) to |. Hence, in the 
following we assume the following simple phenomenological laws ||4^, [47, 4^: 



e = «p^ (23) 

T = (24) 

r = ^— ^ = ap'-\ (25) 

{p + p)cl 

where Cb is the speed of bulk viscous perturbations, i. e., the non-adiabatic contribution to 



the speed of sound in a dissipative fiuid without heat fiux or shear viscosity [|^, and a > 0, 
/3 > 0, r > and s > are constants. Equations (|23|)-(p5D are standard in cosmological 
models, also giving a simple procedure to ensure that the speed of viscous pulses does not 
exceed the speed of light ^ . With this choice we obtain the following scaling equation 
for the bulk viscosity coefficient: 

m = ^ ?(t) = (Lt) , 

Ut) = L'^aLt) = mLt) , (26) 



where we have also used Eqs. (pO|)-([2^). From Eq.(p6D it immediately follows that a scale 
invariant bulk viscosity coefficient must have s = |, and therefore the only physically ac- 
ceptable behavior of the bulk viscosity coefficient is of the form ^ ~ p^/^ ~ . In fact, 
this result immediately follows from Eq. (|26|), since the bulk viscosity coefficient has the 
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same scaling transformation law as the Hubble parameter. In this way, in scale-invariant 
bulk cosmological models the bulk viscosity coefficient is also determined by the expansion 
rate of the Universe. On the other hand, s = | is a case for which the gravitational field 
equations for a flat causal bulk viscous cosmological fluid have been intensively investigated 

0, i9|,|g. 



in the physical literature 



The scaling property of the relaxation time can also be obtained easily and is given by: 



T (t) = ap' 



T (t) 



r (t) 



L^'-^T (Lt) = L- V {Lt) . 



(27) 



The scaling law (p7D holds again only for s = ^. 

In the context of irreversible thermodynamics p, p, T and the particle number density n 
are equilibrium magnitudes which are related by equations of state of the form p = p (T, n) 
and p = p (T, n). From the requirement that the entropy is a state function we obtain 
the equation (|^)^ = ~ ^ {^)n ' "^^^ equations of state (PBD-(P3D, this relation 
imposes the constraint r = (7 — 1) /7 so that < r < l/2forl<7<2, a range 
of values which is usually considered in the physical literature ^7\. Therefore we obtain 
the following scaling law for the temperature of the causal bulk viscous cosmological fluid 
obeying the state equations (|23|-p5D 



T{t) = (3p^ 



{L) 



7— 1 7— 1 



T(t) =(3L'—p— {Lt) 



(L) 

Tit) 



L^"-^T{Lt) 



(28) 



An alternative set of equations of state has been used by Coley, van den Hoogen and 



Maartens [^] to study the full causal dissipative evolution of a FRW space-time. The 

are: 



equations of state proposed in |^ 



1 
H 



X = 



Hx'^, 
P 



(29) 

(30) 
(31) 



with ,^0; and n constants. By writing, with the use of Eqs. (|29| ) - (|3ll) , the Einstein field 
equations as a plane autonomous system, the qualitative behavior and the evolution of the 
Universe can be determined. 
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The scaling behaviors of the bulk viscosity coefficient and of the relaxation time are given, 
for the phenomenological choice of the equation of state Eqs. (p9D-(|3TD, by 

(L) T2m m / r^N 

i it) = S^oLH (Lt) 3^^,4,L (i,) = (Lt) , (32) 

(L) T2n n(T4.\ 

(t) - r ' = LHiLt) ^^^J^l^l^^ = Lr-- {Lt) . (33) 

Therefore the set of equations of state (^9|)-(^) are scale invariant for all m and n. 

For a flat FRW space-time, filled with a perfect fluid obeying a 7-law barotropic equation 
of state, the Hubble parameter H induces a natural time scale via the transformation t — >■ 
in = -^jj- The field equations are invariant with respect to this transformation, having the 
form 3 (^^^)^ = P and 2^^^ + 37i/2 (t^) = 0. Therefore the RG analysis of the 
field equations for the perfect fluid can also be performed in this time variable, and the RG 
transformation in the Hubble time scale tu Ltu leads to equations which are formally 
similar to those obtained by using the usual cosmological time t. The introduction of the 
bulk viscous pressure generally breaks the invariance of the gravitational field equations with 
respect to the Hubble time. In this case, by assuming a reparametrization of the cosmological 
time of the form t ^ to/H, =constant, the 00 component of the field equations is again 
invariant, but the ii, i components take the form (^^^ + ^w) (tn) = 

— H (tn)- This equation is scale invariant only if the bulk viscous pressure is proportional to 
the energy density of the fluid: H ~ = Hop, Hq =constant. Then the numerical coefficient 
in the scaling law of the time is given by to = 2/ (37 + Hq). For all other functional forms of 
the bulk viscous pressure, the Hubble parameter does not generate an invariant time scale 
for the Einstein field equations. 
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IV. RENORMALIZATION GROUP EQUATIONS AND LONG TIME BEHAVIOR 
ANALYSIS 



Letting t = L, from the renormalization group transformations (|n|)-(^) we obtain [|n 

a{t) = t^a(l), 



(t) 



(34) 
(35) 
(36) 

e(t) = r^e(i), (37) 

r(t)=t^r\l). (38) 
To obtain the renormahzation group equations we introduce the parameter of the renor- 



H{t) = r^H{l) 

J*) 

,(*) 



mahzation group transformations in the form L = e^. The infinitesimal transformations of 

(L) (L) (L) 

a , H and 11 , with respect to 6, are given by: 



a a 

~de 

dH 

lie 

dU 



2 (L) da 
2U + 



(39) 
(40) 

de ' ot • (^1) 

With the use of the (scaled) gravitational field equations and of the evolution equation of 
the bulk viscous pressure the renormalization group equations (^ 



d a 

'de 

dH 

'de 

dU 

'de 



2 (L) (L)W 

-— a + a H , 
37 



H + 



1|) can be written as: 
(42) 



2n + 



2 


1{L) 

- -n 
2 




n 

T 


T 


liL) 

- -n 
2 



(i) 1 a™ 

T 



1 



{L)\ 

1 dT 



{L) dt 

T 



J 



(43) 



(44) 



The fixed point (a*, if*, 11*, £*, r*, T*) of the RG equations is given as a solution of the 



system of algebraic equations obtained by setting |11 



da* 



(L) 

dH* 
~df 



(L) 

dU* 
'df 



(45) 
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From the flow in the RG around the fixed point we can see the long time-behavior of 
causal bulk-viscous fiuid filled fiat Universe, once the equations of state for the bulk viscosity 
coefficient, relaxation time and temperature are given. 

As an application of the mathematical formalism presented above we shall consider the 
long time-behavior of the causal bulk viscous cosmological model obtained by adopting the 
equations of state (P^)-(PT|), with s = | and r = In this case the gravitational and 

causal bulk evolution equations are scale-invariant. The equations of state take the simple 
form ^ = \/3aH, r = and T = The causal bulk viscous pressure evolution 

equation is 



rn + n 



-3^H - -tU 



3if-2(l + r)| 



(46) 



while the RG equation (|4^) becomes 



dU 



(L) 

2n + 



n p (L) liL) / (L) n 

__3X//-5n 6,ff + (l + r)- 

^ ^ \ H. 



(47) 



The fixed point of the RG equations 



and ( ^7\ ) is defined by 



4 J ' 37'^^ 372' 



0,r = 0,r 



0. 



(48) 
(49) 



The fixed point of the RG equations corresponds to a fiat, non- viscous Friedmann Uni- 
verse. In order to study the fiow in the RG around the fixed point, we must study the 

(i) (L) (L) (L) 

perturbation around the fixed point. We define the perturbed quantities 5 a , 5H , 511, 6 ^ 
and (5r-i by 111 



a 



a + a , 



(L) 

H =H' 



6H, 



n = u* + 6U, 

(L) (L) (L) 

^ =^* + S^ =V3aSH, 

(L) {L) /o (L) 

a 



(50) 
(51) 
(52) 
(53) 
(54) 
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(i) (L) (L) (L) (i) 

where the perturbations 5a, (5if, 511, 5 ^ and 5r are assumed to be small quantities. 



6 a « 1 etc. From the RG equations (|42|), (^31) and p7|), it follows that the perturbed 



quantities satisfy the following linearized differential equations 

a*6H, (55) 



d9 

dSH W 1 W 

^^-5i/--5n, (56) 

From Eqs. (|55| ) and ( ^61 ) we obtain 

while Eqs. ( |5B| ) and ( |57D give the following evolution equation for the perturbed bulk viscous 
pressure: 

(L) [L) 

d'^5Tl d5Tl 2/1 1\ W , , 

+ -^--H^--Mn=0. (59) 



dO'^ de 7 VTSa 7 
The general solution of Eq. (ESI) is given by 



5Ii {6) = + C2e''^\ (60) 



with Ci, C2 constants of integration and 



-•i,! = ^- 1^ (61) 

From (|6TD it follows that r2 < 0, Va, 7 while the sign of ri depends on the sign of the 
quantity ~ 7) ' ^'^ (^fe "7) > 0, ( 7 > v^a ) then ri > while for 7 < v^a, ri < 0. 

With the use of Eq. (|60|) , equation (^), describing the evolution of the perturbation in 
RG of the scale factor of the causal bulk viscous fluid filled Universe takes the form 

2 (^) (^) 

and has the general solution given by 

5^? {d) = he-' + he^^' + he^-\ (63) 
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where fi,i = 1,2,3 are constants of integration. 

From Eq. ( |63D we can see the flow in the RG around the fixed point. A fixed point yo 
of a system y = F{y) is said to be stable if for every neighbourhood U of yo, there is a 
neighbourhood U' < U oi y^ such that every trajectory which passes through U' remains 
in f/ as t increases. The fixed point is asymptotically stable if it is stable and there is a 
neighbourhood U of yo such that every trajectory passing through U approaches y^ as t 
tends to infinity. A fixed point which is not stable is said to be unstable ||53| . 

For 7 < -y/Sa the fixed point is an attractor. Then, by setting the initial profile in the 



vicinity of the fixed point (^81) - (^9|), the fiat causal bulk viscous Universe with bulk viscosity 
coefficient proportional to the Hubble parameter will approach, in the long-time limit, the 
fiat Friedmann Universe with a{t) = a*t^/'^^ . This point is also asymptotically stable. For 
7 > v^tt, there is a single relevant mode. Therefore, if 7 > -y/Sa, the space-time will deviate 
from the fiat Friedmann geometry since there is a relevant mode 5a{t) = f^t^^- This mode 
corresponds to a growing mode, and therefore in this case the fixed point is asymptotically 
unstable. 

The conditions of the asymptotic stability or instability of the solution can also be for- 
mulated in terms of the speed of sound in the cosmological fiuid. Since \/3a = ^/H = 
VS^/y/p = v^3p7c^r, the condition of asymptotic stability becomes > l/^/3pT. With the 
use of the relation r = ap~^^'^, the conditions of asymptotic stability or instability are given 
by cl > l/VSa and cl < I/a/So, respectively. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have performed an analysis of the fiat homogeneous bulk viscous 
cosmological models by using the RG method, in the version developed for the study of 
asymptotics of non-linear evolution equation. For the fiat FRW geometry the Einstein grav- 
itational field equations are invariant with respect to the scale transformations of the time, 
energy density and Hubble parameter. Therefore it is required that the causal bulk evolution 
equation must also be scale invariant with respect to the same group of transformations,thus 
leaving invariant the full set of gravitational field and bulk viscous pressure equations. This 
condition leads to specific scale transformation laws for the bulk viscosity coefficient and 
relaxation time, and imposes strong constraints on the equations of state relating these pa- 



16 



rameters to the energy density of the cosmological fluid. By using the property of the scale 
invariance of the field equations, we have shown that the most used system of equations of 
state, expressing the bulk viscosity parameter and relaxation time as a power of the energy 
density, ~ p'*, r ~ p*^^, is scale invariant for only one value of the exponent s. Of course 
this result crucially depends on the linear barotropic equation of state adopted to describe 
the matter content of the very early Universe. This equation of state is not satisfied in more 
realistic physical cases, in which, for example, the Maxwell-Boltzmann gas could provide a 



more appropriate description of the cosmological fiuid ||25[| . 

The RG analysis also provides a powerful method for the study of the long-time asymp- 
totics of the gravitational field equations. This can be done by the study of the fiow around 
the fixed point of the RG group equations. For the cosmological model having the bulk 
viscosity coefficient proportional to the Hubble parameter H, the fixed point corresponds to 
a fiat ideal (non-viscous) Friedmann type geometry. Due to the presence of viscous effects 
the geometry of the perturbed Universe will deviate from that of the fixed point, and the 
linear analysis of perturbations indicates the final state of the bulk viscous Universe. In the 
considered cosmological model the long-time behavior of the space- time depends on the 
ratio of the constant a relating ^ to H and 7. 

Bulk viscous cosmological models have been investigated in the past years from many 
different points of view. The method of the renormalization group, as applied in a systematic 
manner to this class of models, can become a valuable tool for cosmologists leading to a better 
understanding of the nature, behavior and properties of our Universe. 
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